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Abstract 

The evolution of a modulated positron beam in a planar crystal channel 
is investigated within the diffusion approach. A detailed description of the 
formalism is given. A new parameter, the demodulation length, is introduced, 
representing the quantitative measure of the depth at which the channelling 
beam preserves its modulation in the crystal. It is demonstrated that there 
exist crystal channels with the demodulation length sufficiently large for using 
the crystalline undulator as a coherent source of hard X rays. This finding is 
a crucial milestone in developing a new type of lasers radiating in the hard 
X ray and gamma ray range. 
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1. Introduction 

In this article we study the evolution of a modulated positron beam in 
straight and bent planar crystal channels. Some key ideas of this research 
were briefly communicated in [l[ and jij . In this paper we present a system- 
atic and detailed description of the formalism and the obtained results. The 
outcome of the research is of crucial importance for the theory of the crystal 
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undulator based laser (CUL) 0, 0, Isf — a new electromagnetic radiation 
source in hard x- and gamma-ray range. 

Channelling takes place if charged particles enter a single crystal at small 
angle with respect to crystallographic planes or axes joj. The particles get 
confined by the interplanar or axial potential and follow the shape of the 
corresponding planes and axes. This suggested the idea of using bent 
crystals to steer the particles beams. Since its first experimental verification 
jsj the idea to defiect or extract high-energy charged particle beams by means 
of tiny bent crystals replacing huge dipole magnets has been attracting a 
lot of interest worldwide. Bent crystal have been routinely used for beam 
extraction in the Institute for High Energy Physics, Russia A series of 
experiments on the bent crystal defiection of proton and heavy ion beams was 



performed at different accelerators [10|, llll, |12|, ll3|, [l^ throughout the world. 



The bent crystal method has been proposed to extract particles from the 
beam halo at CERN Large Hadron Collider 15| The possibility of defiecting 
positron 16|] and electron 



14, 17 beams has been studied as well. 




Figure 1: Schematic representation of the crystahine undulator. 

A single crystal with periodically bent crystallographic planes can force 
channelling particles to move along nearly sinusoidal trajectories and radiate 
in the hard x- and gamma-ray frequency range (see figure [1]). The feasibility 



2 



of such a device, known as the 'crystalhne undulator', was demonstrated 
theoretically a decade ap;o (further developments as well as historical 



references are reviewed in |18|). M ore recently, an electron based crystalline 
undulator has been proposed [l9|]. 

It was initially suggested to obtain sinusoidal bending by the propagation 
of an acoustic wave along the crystal js], 0] • The advantage of this approach 
is its flexibility: the period of deformation can be chosen by tuning the fre- 
quency of the ultrasound. However, this approach is rather challenging tech- 
nologically and yet to be tested experimentally. Several other technologies 
for the manufacturing of periodically bent crystals have been developed and 
tested. These include making regularly spaced grooves on the crystal surface 
either by a diamond blade 20, 21 1 or by means of laser- ablation 22[, depo- 
sition of periodic Si3N4 layers onto the surface of a Si crystal 21 1 , g rowing 
of Sii-i^Gci: crystals[23] with a periodically varying Ge content x |2J, |25 . 

Experimental studies of the crystalline undulator are currently in progress. 
The first results are reported in 



26 and 27 



The advantage of the crystalline undulator is in extremely strong electro- 
static fields inside a crystal which are able to steer the particles much more 
effectively than even the most advanced superconductive magnets. This fact 
allows to make the period Au of the crystalline undulator in the range of hun- 
dreds or tens micron which is two to three orders of magnitude smaller than 
that of conventional undulator. Therefore the wavelength of the produced 
radiation A ~ Au/(27^) (7 ~ lO'^-lO^ being the Lorentz factor of the par- 
ticle) can reach the (sub)picometer range, where conventional sources with 



comparable intensity are unavailable |28 



Even more powerful and coherent radiation will be emitted if the prob- 
ability density of the particles in the beam is modulated in the longitudinal 
direction with the period A, equal to the wavelength of the emitted radia- 
tion (see figure |2]). In this case, the electromagnetic waves emitted in the 
forward direction by different particles have approximately the same phase 



Therefore, the intensity of the radiation becomes proportional to the 
beam density squared (in contrast to the linear proportionality for an un- 
modulated beam). This increases the photon flux by orders of magnitude 
relative to the radiation of unmodulated beam of the same density. The ra- 
diation of a modulated beam in an undulator is a keystone of the physics of 
free-electron lasers (PEL) (sT 



31 



It can be considered classical coun- 
terpart of the stimulated emission in quantum physics. Therefore, if similar 
phenomenon takes place in a crystalline undulator, it can be referred to as 
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Figure 2: In an unmodulated beam (the upper panel) the particles are randomly dis- 
tributed. In a completely modulated beam (the lower panel) the distance between any 
two particles along the beam direction is an integer multiple of the modulation period A. 



the lasing regime of the crystalline undulator. 

The feasibihty of CUL radiating_in the hard x-ray and gamma-ray range 
was considered for the fist time in |3|, 4]. Recently, a two-crystal scheme, the 
gamma klystron, has been proposed 5|]. 

A simplified model used in the cited papers assumed that all particle tra- 
jectories follow exactly the shape of the bent channel. In reality, however, the 
particle moving along the channel also oscillates in the transverse direction 
with respect to the channel axis (see the shape of the trajectory in figure 
[1]). Different particles have different amplitudes of the oscillations inside the 
channel (figure [3l upper panel). Similarly, the directions of particle momenta 
in [xz) plane are slightly different (figure El lower panel). Even if the speed 
of the particles along their trajectories is the same, the particles oscillat- 
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Figure 3: Due to different amplitudes of channelling oscillation (upper panel) and different 
momentum directions in the (xz) plane (lower panel), the initially modulated beam gets 
demodulated. The open and filled circles denote the same particles at the crystal entrance 
and after travelling some distance in the crystal channel, respectively. 

ing with different amplitudes or tlie particles with different trajectory slopes 
with respect to z axis have slightly different components of their velocities 
along the channel. As a result, the beam gets demodulated. An additional 
contribution to the beam demodulation comes from incoherent collisions of 
the channelling particles with the crystal constituents. 

In the case of an unmodulated beam, the length of the crystalline undu- 
lator and, consequently, the maximum accessible intensity of the radiation 
are limited by the dechannelling process. The channelling particle gradu- 
ally gains the energy of transverse oscillation due to collisions with crystal 
constituents. At some point this energy exceeds the maximum value of the 
interplanar potential and the particle leaves the channel. The average pen- 
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etration length at which this happens is known as the dechannelling length. 
The dechannelled particle no longer follows the sinusoidal shape of the chan- 
nel and, therefore, does not contribute to the undulator radiation. Hence, 
the reasonable length of the crystalline undulator is limited to a few dechan- 
nelling lengths. A longer crystal would attenuate rather then produce the 
radiation. Since the intensity of the undulator radiation is proportional to 
the undulator length squared, the dechannelling length and the attenuation 
length are the main restricting factors that have to be taken into account 
when the radiation output is calculated. 

In contrast, not only the shape of the trajectory but also the particles 
positions with respect to each other along z axis are important for the lasing 
regime. If these positions become random because of the beam demodulation, 
the intensity of the radiation drops even if the particles are still in the chan- 
nelling mode. Hence, it is the beam demodulation rather than dechannelling 
that restricts the intensity of the radiation of CUL. Understanding this pro- 
cess and estimating the characteristic length at which this phenomenon takes 
place is, therefore, a cornerstone of the theory of this new radiation source. 

2. Diffusion Equation 

2.1. The model of the crystal channel 

We adopt the following model of the planar crystal channel (see Fig. 12. ip : 
- the interplanar potential is approximated with a parabola 



(p is the distance from the potential minimum) so that the channeling oscil- 
lations are assumed to be harmonic; 

- the electron density within the distance of one Thomas-Fermi radius of the 
crystal atoms from the crystallographic plane is assumed to be so high that 
the particle gets quickly scattered out of the channel. Therefore, the particle 
is considered dechanneled just after it enters this region. So that the effec- 
tive channel width is 2pmax = d — 2aTF, where d and otf are respectively the 
interplanar distance and is the Thomas-Fermi radius. 

As is seen from the figure, the parabolic approximation is quite reason- 
able. The real potential differs from the parabola mostly in the region of 
high electron density, where the particle assumes to be dechanneled. 




(1) 
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Figure 4: The model of the plane crystal channel. The interplanar potential is approxi- 
mated by a parabola. It is assumed that the particle dechannels if it enters the vicinity of 
the crystallographic plane within the Thomas- Fermi radius, atf- 



2.2. The particle distribution 

Let us consider the distribution f{t,s;^,Ey) of the beam particles with 
respect to the angle between the particle trajectory and axis z in the (xz) 
plane ^ = arcsmpx/p ~ Px/p and the energy of the channeling oscillation 
Ey = Py/2E + U{y) 0. Here p, p^ and Py are, respectively, the particle 
momentum and its x and y components, and E is the particle energy (we 
will consider only ultrarelativistic particles, therefore E ^p). 



^We chose the system of units in such a way that the speed of light is equal to unity. 
Therefore, mass, energy and momentum have the same dimensionality. This is also true 
for length and time. 
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2.3. Kinetic equation 

In absence of random scattering, the distribution function f{t,z;^,Ey) 
would satisfy the differential equation ^ + = 0, where Vz = ||. In 
reality, however, the right-hand-side of the equation is not zero. It contains 
the collision integral. After averaging over the period of the channeling os- 
cillation, the kinetic equation takes the form 

% + %(^^) = Qjd^dE'y[f{t,z;e,E'^)w{e,E'^;^,Ey) 

-f{t,z;tEy)w{^,Ey;e:K)]) (2) 

where w{^,Ey;^',Ey)dz is the probability that the particle changes its angle 
^ and transverse energy from ^ and Ey to, respectively, ^' and Ey while 
travelling the distance dz. The angular brackets stand for averaging over the 
period of the channeling oscillations. 
Due to the detailed equilibrium 

w{tEy;e:E'y)^w{e,E'y,^,Ey) (3) 

2.4- Diffusion approximMion 

We assume that soft scattering dominates, i.e. the function to (^', Ey \ ^, Ey) 
is not negligible only if — ^| and \Ey — Ey\ are small so that /(t, z] Ey) = 
f{t, z;^ + dx, Ey + Qy) can be expanded into the Taylor series with respect to 
-dx and Qy. Then, up to the second order in -d^ and Qy, one obtains 

df df , , ^df ^ df ^ d^f ^ d^f ^ d'^f 



where 



(5) 



-D^ = (^J d'dx'dx J dqyW{C,Ey;^ + -dx,Ey + qy) 

Dy = l^j d'dx J dqyqyW{E^,Ey;^ + i&x,Ey + qy)'^ (6) 

= IQ € J dqyw{i, Ey- ^ + i?,, Ey + g,)^ (7) 

D^y = (^j d'dx'dx j dqyqyw{^,Ey]i + 'dx,Ey + qy)^ (8) 

Dyy ddx j dqy q^wiC, Ey-, ^ + ^x, Ey + g^)^ (9) 
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3. Diffusion Coefficient 



3.1. Scattering 

Let us consider a channeling positron colliding with a target electron. If 
6 is the scattering angle in the lab frame and </? is the angle between the 
scattering plane and the (a;2;)-plane then the transverse components of the 
particle momentum are changed by 

^Px — psm9cos(p, (10) 
^Py — p sin ^ sine/?, (11) 

As far as 6* <^ 1, we can use the approximation sin 6' ~ 9. Then 

q Px + Spx Px hx a no^ 
Vx — = — = 9 COS (p. (12) 

P P P 



and 



Qy = 



Py9smip + -9 sm (p. 



3.2. The transition probability 

The probability for the particle to be scattered by an electron from the 
state (^, Ey) to the state (^ + '&x, Ey + Qy) while travelling the distance dz can 
be related to the differential cross section of positron-electron scattering: 

w{^,Ey]^ + ^x,Ey + qy)dz = nedz J d9 J ^(p-^^ (14) 

5 {9 cos ip — '&x) S (py9 sin ip + ^9"^ sin^ cp — qy^ 

Because both target and projectile are not polarized, the cross section does 
not depend on ip: 

^ = J-^. (15) 
d9d(p 2Tr d9 ^ ^ 
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Substituting (fn|) into ([S])-© and integrating over -d^ and qy one obtains 



= ^y'ejdej^ dif^^pyesmif + ^e^sm^ip)) (17) 

1 / /" /"^"^ tic'" \ 



Df„ = — (fie / c?^'^ / (iv96' COS (y9 (]9„6' sin H — 6'^sin^y9)) (19) 
2-n \ J dO Jq \ 2 / / 

= ^(^-/^^^/ rf¥^^^inV(p, + |^sin¥p)'^ (20) 

Then integration over and neglecting higher order terms with respect 
to 9 yields 

= (21) 

Dy = I de^e' (22) 

D,, = \{n.) I de^e' (23) 

= (24) 



1 , o, f ,„da 



Dyy = -^{UeVD J ^^^^ (25) 

Here (ue) is the electron density along the particle trajectory averaged 
over the period of the channeling oscillations. Generally speaking, {rie) de- 
pends on the transverse energy Ey We assume, however, that the electron 
density does not change essentially within the channel. Therefore, (ng) can 
be treated as a constant. For the same reason, we can make the approxi- 
mation (uepl) ~ {n'e){pl) Then {py) = 2E (^^^ = EEy. due to the virial 

theorem for the harmonic potential: Py/(2E) = Ey/2 
Finally, one obtains for nonzero coefficients 



Dy = Do (26) 
^« = (27) 

Dyy = EyDo (28) 
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The diffusion equation talces the form 



df ^df 



d 



BE,. 



E, 



dl 
'dE„ 



+ 



EOe 



(29) 



Equation fl2^ is akin to the equation describing dechannehng process (see 
e.g. I32I). The novel feature of it is the presence of time variable, which 
allows to describe time dependent (modulated) beams. Additionally, it takes 
into account scattering in the {x, z) plane. 



4. Solving the diffusion equation 

4-1. The average longitudinal velocity 

The particle velocity along z axes averaged over the period of channeling 
oscillations can be represented as 



'1 - ^ cos i 

(v.) = ' (30) 

Here a/1 — 1/7^ is the particle speed along the trajectory, cos^ ^ (1 — ^^/2) 
appears because of the slope ^ <^ 1 of the trajectory to z axis in (xz) plane, 
and the denominator is due to the sinusoidal channeling oscillations in (xy) 
plane with the amplitude b and the period Ac = 2TT/kc. Taking into account 
that the amplitude of the channeling oscillations is much smaller than their 
period, bk^ ^ 1, the denominator can be approximated by 1 + {bk^Y /A. For 
the harmonic potential ([I]) (see Fig. 12.11) the amplitude b is related to the 
transverse energy Ey by 

& = PmaxyT^. (31) 

Using the formula for the frequency of the harmonic oscillator one finds 



So that bkc = \l • Finally, neglecting higher order terms 

{v.) ^ ( 1 - At - - - ^1 (33) 
^ ' ^ 272 2 2Ey ^ ^ 
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4-2. Excluding the time variable 

If the beam is periodically modulated (bunched) the distribution f{t, z; ^, Ey) 
can be represented Fourier series: 



f{t,z;^,Ey)= J2 

9j{z;^,Ey)exp{ijut). 



(34) 



j=-oo 



with g*{z; ^, Ey) = g-j{z; ^, Ey) to ensure the real value of the particle distri- 
bution. Since Eq. fl29|) is linear, it is sufficient to consider only one harmonic. 
Substituting f{t, z; ^, Ey) = g{z; ^, Ey) exp^iut) into f l29l) one obtains 



dg 

i^9{z] Ey) + = Do 



d 



dE„ 



E,, 



dg 
BE,, 



1^ 
Edi^ 



(35) 



4-3. Variable separation 

To simplify this equation, we make the substitution 

^(^; Ey) = exp {-iujz) g{z; ^, Ey), 

where g{s; ^, Ey) varies slowly comparing to exp {—iuz): 

dg/dz < ujg{z;^,Ey). 

Equation ( l35l) takes the form 



dg_ 

dz 



Vz) + iojg{z] Ey){l - {v^)) = Do 



d f ^ dg 



dEy V "^dEy 



+ 



1^ 

Edi^ 



(36) 



(37) 



(38) 



In the first term, the velocity can be approximated by unity: (vz) ~ 1, i.e. the 
term dg/dz{l — (vz)) can be neglected. However the term iug{z; C,, Ey){l — 
(vz)) has to be kept because of (1371) . Using the expression (l33l) for (vz), one 
obtains from f l38|) the following partial differential equation for g{z] ^, Ey) 



dg{z]^,E, 



dz 



lUJ 



d 



" +i7^9{z;i,Ey) = Do^\E 

K.^, . Dod^~g{z;^,Ey) e 



•dEy^-' 



dg{z; ^, Ey 



dE„ 



(39) 



—i(^ — Qiz;C,E^.)+ ^ ^ 



^^i^gi.z]i,Ey) 
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This equation can be solved by the method of separation of variables. Putting 
g{z;$,,Ey) = Z{z)E{^)£{Ey), after substitution into we obtain a set of 
ordinary differential equations: 

Do 1 d'Ejo . e . .... 

= (40) 



E{Ey) dEy \ ' dEy J 2E 
1 dZ{z) ico 
Z{z) dz ^2^ ^ ^ ' 

where C^, Q and Cy do not depend on any of the variables z, ^ and Ey and 
satisfy the condition 

C, = C^ + Cy. (43) 

4-4- Solving the equation for 

Equation fHOj) can be rewritten as 

We change the variable 



and introduce the notation 



This results into 



^^-X'^{x) = -mx). (47) 

This equation has the form of the Schrodinger equation for the harmonic 
oscillator. Its eigenvalues and integrable eigenfunctions are well known: 

Qn = 2n + l (48) 
S„(x) = i/„(x)exp(-xV2), (49) 
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where n = 0, 1, 2, . . . and Hn{x) = ("f^) ^ Hermite Polynomials 

satisfying the orthogonality condition 

/ + 00 
dx e-^"Hn{x)Hn'{x) = 5nn'2"n!0F (50) 
-oo 

which is equivalent to 

/ + 00 
dxE4x) 2„'(X) = 5nn'2"n!v^. (51) 
-oo 

Returning back to the variable ^ one obtains 



H„«).^„|e-y^,je.p(-i±iyff|, (52) 

Any integrable function F{C,) can be represented as series 



FiO = F [e-^'f^xj = E ^n^r^ix) (53) 
Let us multiply the above expression by S„/(x) and integrate over x 



/+00 / Ion \ °° r+co 



iO (54) 



Using fl5U]) one finds 



2"n!v^ 
From (HHI) and fHHl) one finds 



^ / °°rfxe^-^/^F e--/^y^x ^^n(x) (55) 



Q,„ = -(l + z)'^/^U + i), n = 0,l,2,... (56) 
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4-5. Solving the equation for £{Ey) 
Equation fHT]) can be rewritten as 



' ^E;^)J^E,,^]aE,)-0, (57) 



dEy \ ' dEy J \2DoE ' Dq) 
By the substitution 



equation (l57Il can be reduced to 



Ey = ^—J^e (58) 



a; 



with 



2.+ l = -(l->)y-|-C, (60) 

Further substitution £{e) = exp{—6/2)L{e) results into the Laguerre 
equation: 

One of two hnearly independent solutions of this equation is logarithmi- 
cally divergent at e — )■ and, therefore, has to be rejected. Another solution, 
Ly{e)^ is finite at e = and is known as the Laguerre function 

Returning back to the variable Ey, the solution of equation (jS7|l can be 
represented as 

£(B.)=exp(-i±iy£:^„)^(l + ,)y^B,) (62) 

The eigenvalues can be found by imposing the boundary conditions. If 
the energy of the channeling oscillations exceeds the value ?7max (see Fig. 12. ip 
the particle enters the region of high electron density, get scattered by crystal 
constituents and becomes dechanneled. Therefore, the distribution function 



^ At nonnegative integer values of v, the Laguerre function is reduced to the weh known 
Laguerre polynomials. In the general case that is relevant to our consideration, it can be 
represented by an infinite series ( [Appendix A. 12 1. 
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of channeling particles has to be zero at Ey = ?7max- This results into the 
following boundary condition 



Equation fl63l) has to be solved for u. Then, according to ( 160|) . the eigenvalue 
Cy^k can be found from 



[l + t) iDou^ 
E 



Cy = -^^\l^i2u+l). (64) 



The subscript k = 1, 2, 3, . . . enumerates different roots of equation (163|1 . 
We introduce a dimensionless parameter 

^ UL. (65) 



jiiDoE 

{jo,k is k-th zero of the 0-th order Bessel function: Jo(jo,fc) = 0). Then 
equation flB51) can be rewritten as 

^Jo.iv^l = 0. (66) 



This equation has infinite number of complex roots (see Appendix) which we 
denote as //^(k), = 1, 2, 3, . . . . The equation does not have any analytical 
solution and therefore has to be solved numerically. 

Instead of the complex function //^(k), it is more convenient to introduce 
two real functions: 



a.in) = ^[l + 2muk{K)]-Q[ukiK)])] (67) 
Jo,i 

Pk{^^) = -^[l + 2(3fJ[z/,(/«)] + 5>[z/fe(K)])]. (68) 
^Jo,iv 

The eigenvalues can be represented in the form 

Cy,k = -^-^u:9l^3,{K). (69) 

Here 

Ld = ^Um^JUlM and (70) 
= ^2U^.JE (71) 
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are, respectively, the dechanneling length 32|] and Lindhard's angle. The 
parameter k fl65l) can be rewritten in terms of and 



where A = 27r/u; is the spatial period of the modulation. 

4-6. Solving the equation for Z{z) 
Equation has the solution 



(72) 



Z{z) = exp ( C, - i— 



(73) 



The value of Cz can be found using (H3l) . ( 1561) and (JMI) . Then the solution 
(1731) takes the form 



exp <; - — 



afc(K) + (2n + 1)^ 
Jo,i 



(74) 



lUJZ 



1 



27 2jo,iv^ 



Hence, the solution of Eq. ( I35|) is represented as 

oo oo 

5f(2:; ^, ^j^) = exp {-iuz) ^ ^ an,k'^n{0^k{Ey)Zn,k{z), 



(75) 



n=0 fc=l 



where the coefficients an,k are found from the particle distribution at the 
entrance of the crystal channel: 



^n,k 



(76) 



5. The demodulation length 

5.1. The demodulation length in a straight channel 

Due to the exponential decrease of Zn,k{z) with z (see (17^ ). the asymp- 
totic behaviour of g{z; ^, Ey) at large z is dominated by the term with n = 
and k = 1 having the smallest value of the factor [^^(k) + {2n + l)y^/jo,i] 
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in the exponential. Therefore, at sufficiently large penetration depths, the 
particle distribution depends on z as 



g{z\i,Ey) ocexp(-z/Ldm -iuj/u^z) (77) 
where Ldm is the newly introduced parameter — the demodulation length: 



^dm 



(78) 



and Uz is the phase velocity of the modulated beam along the crystal channel 



2 ' 



27 V 2jo,iVK 



(79) 



This parameter is important for establishing the resonance conditions be- 
tween the undulator parameters and the radiation wavelength. 

In this article we concentrate our attention on the demodulation length. 
This parameter represents the characteristic scale of the penetration depth 
at which a beam of channeling particles looses its modulation. 

Fig. 15. II presents the dependence of the ratio Ldm/-^d on the parameter k. 
It is seen that the demodulation length approaches the dechanneling length 
at K < 1. On the contrary, the ratio noticeably drops for k > 10. 

It is instructive to study the influence of the particle motion in x and 
y direction on the demodulation length separately. Replacing ai{K) in ( !78l) 
with unity means neglecting the motion in the y direction, while omitting the 
second term in the denominator ignores the motion in x direction. One sees 
from Fig. 15. II that it is mostly the motion in x direction that diminishes the 
demodulation length at k < 10, while the influence of channeling oscillations 
is negligible. This suggests the idea that for the axial channeling, i.e. when 
motion in both x and y directions has the nature of channeling oscillations, 
the demodulation length L^m may practically coincide with the dechanneling 
length Ld at higher frequencies of the beam modulation than in the case of 
planar channeling. 

5.2. The centrifugal force in a bent channel 

So far, beam demodulation in a straight channel has been considered. 
The channels of a crystalline undulator, however, have to be periodically 
bent. Therefore the above formalism has to be modified for the case of a 
bent channel. 
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Let us consider a crystal that is bent in the (yz) plane so that the crystal 
channel has a constant curvature with the radius R. An ultrarelativistic 
particle with energy E moving in such a channel experiences the action the 
centrifugal force 

^c.f. = |. (80) 

It is convenient to characterise the channel curvature by the dimensionless 
parameter C defined as 



C 



U' 

max 



(81) 



where U'^^^ is the maximum value of the derivative of the particle potential 
energy in the channel, i.e. the maximum transverse force acting on the 
particle in the interplanar potential. Channeling is possible at < C < 1. 
The value C = corresponds to a straight channel. The critical radius Rc 
(known also as Tsyganov radius) at which the interplanar potential becomes 
unable to overcome the centrifugal force corresponds to C = 1. 
In the case of potential energy ([T]), 

UL. = f/'(Pmax) = 2^. (82) 
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so that 



C 



Pmax-^c.f . 

OTJ 

max 



(83) 
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Figure 6: The potential energy of a particle in the planar crystal channel for a straight, 
C = 0, and for a bent, C 7^ 0, crystal. The effective width of the bent channel is 
2pmax(l — C) and the depth of the potential well is L/max(l — C)'^, where 2pmax and f/max 
are, respectively, the effective width and the depth of the straight channel (cf. Fig. 12. ip . 



The potential energy is modified by the centrifugal force in the following 

way 

Uc{p) = U{p)~pF,,,. (84) 
For the parabolic potential energy ([1]) the modified potential can be conve- 



20 



niently rewritten in terms of the parameter C\ 



Uc{p) = f/n 



2 

2 



C\ -C 



(85) 



The potential energy Uc{p) reaches its minimum at po = Cpmax- The effective 
width of the channel becomes (see Fig. 15. 2p 

Pmax - Po = Pmax(l " C). (86) 

The depth of the potential energy well is 

f/c(Pmax) - Ucipo) = f/max(l - Cf. (87) 

So to obtain the solution of the diffusion equation for the bent crystal we can 
use the results of Sec. H] with the substitution 

f/max^f/max(l-C)2. (88) 

5.3. The demodulation length in a bent channel 

Substitution (l88ll modifies the demodulation length and the Lindhard's 
angle the parameter k in the following way: 

Ld -> L^l-C)^ (89) 
^ e^{l-C) (90) 

Consequently, the the modification of parameter k is 

k^k(I-C)^ (91) 

It is convenient to introduce modified functions afc(/t, C) and f3k{K,C): 

Mf^^^) = — _ — (92) 
/3fc(/€,C) = (l-C)^/3,.(/€(l-C)^). (93) 
In terms of these functions, the eigenvalue Cy^k has the form 

Cy,k = -^^^^^ - zujel/3k{K, C). (94) 
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This exactly coincides with ( 1^ up to replacing akin) and /^^(k) with a;fc(/t, C) 
and Pk{i^, C), respectively. Note that and 6'l in fl94p have the same mean- 
ing as in fl69|) : they are related to the straight channel. 

Similarly, the demodulation length in the bent channel is given by 



^dm 



(95) 
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Figure 7: The ratio of the demodulation length Ldm ([95)) to the dechannelhng length in the 
straight channel Ld versus the parameter k (|72p for different values of curvature parameter 
C. The corresponding asymptotic values at k ^ are shown by thin horizontal lines. 



Fig. 15.11 presents the dependence of the ratio Ldm/-^d on the parameter 
K. At K — )■ 0, the demodulation length approaches (1 — C^Ld which is 
the dechannelhng length in a bent crystal. It is seen that the demodulation 
length is smaller than dechannelhng length by only 20-30% at k < 1 for C 
ranging from to 0.3. It noticeably drops, however, at k > 10. 

The above estimations are made for the channel with constant curvature. 
Similar results are expected for a sinusoidal channel with C varying between 
and 0.3. 

It was proven for a number crystals channels [HI that the dechannelhng 
length of positrons is sufficiently large to make the crystalline undulator 
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feasible. Such a crystalline undulator becomes a CUL, i.e. it generates 
coherent radiation, provided that it is fed by a modulated positron beam 
and the beam preserves its modulation over the length of the crystal. This 
takes place if the demodulation length in the crystalline undulator is not 
much smaller than the dechannelling length. Hence, CUL is feasible if there 
exist crystal channels ensuring k < 1 in the range of the photon energies 
above ~ 100 keV (softer photons are strongly absorbed in the crystal). It 
will be shown in the next section that such crystal channels do exist. 

5.4- Estimation of the parameter k 

To evaluate the parameter k (172|) we shall use the approximate formula 



for the dechannelling length |32| : 



256 E otf d 
^'^ = in T- ^^^> 

Here and tq are, respectively, the electron mass and the classical radius, 
d is the distance between the crystal planes, and the Coulomb logarithm A 



for positron projectiles is defined as [18 



A = ,og^H^-5, (97) 

& J 24' ^ ' 

with 

/ ^ 16Z°-^ eV (98) 

being the ionization potential of the crystal atom with the atomic number 
Z. The Thomas- Fermi radius of this atom is related to the Bohr radius ob 
by the formula 

0-8853 

otf = Qb 3^ • (99) 
V Z 

Substituting (jnSD into (1721) and taking into account fITT]) . one obtains 

512 OtF f/max d 



Qtt Aro me A 



(100) 



As is seen from the above formula the value of k is determined by the potential 
depth t/max, by the distance between the planes d and the modulation period 
A. It also depends on the atomic number of the crystal atoms Z via (1981) and 
These parameters are listed in table 15.41 for several crystal channels. 
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Table 1: The parameters of the crystalchannels used in the calculations (see text). For 
(111) plane of diamond, only the larger of two channels is presented. 



Crystal 


Z 


/(eV) 


axF (A) 


Plane 


d (A) 


f/max (eV) 


Diamond 


6 


80 


0.26 


(100) 


0.9 


2.2 










(110) 


1.3 


7.3 










(lll)L 


1.5 


10.8 


Graphite 


6 


80 


0.26 


(0002) 


3.4 


37.9 


Silicon 


14 


172 


0.19 


(100) 


1.4 


6.6 










(110) 


1.9 


13.5 


Germanium 


32 


362 


0.15 


(100) 


1.4 


14.9 


Tungsten 


74 


770 


0.11 


(100) 


1.6 


56.3 



The dependence on the particle energy is cancelled out, except the weak 
dependence due to the logarithmic expression (p7I) . 

The dependence of the parameter k on the energy of the emitted photons, 
hjj = 2TTh/X, is shown in Fig. 15.41 The calculation was done for 1 GeV 
positrons. Due to the weak (logarithmic) dependence of k on the particle 
energy, changing this energy by an order of magnitude would leave Fig. 15.41 
practically unaltered. 

As one sees from the figure, k ~ 1 corresponds to /jo; = 100 — 300 keV 
for (100) and (110) planes in Diamond and (100) plane in Silicon. So these 
channels are the most suitable candidates for using in CUL. This is, however, 
not the case for a number of other crystals e.g. for graphite and tungsten 
having k > 10 in the same photon energy range. 

At /zo; ~ 10 MeV, k becomes larger than 10 for all crystal channels. This 
puts the upper limit on the energies of the photons that can be generated by 
CUL. It is expected to be most successful in the hundred keV range, while 
generating MeV photons looks more challenging. 

6. Discussion and Conclusion 

One may expect that the demodulation is not limited to the processes 
illustrated in figure [31 An additional contribution can come from the energy 
spread of the channelling particles, as it usually happens in ordinary FELs. 
In fact, the contribution of the energy spread to the beam demodulation 



24 



10000 



1000 



100 



10 



0.1 



0.01 



C(diam.)(100) 
C (diam.) (110) 
C (diam.) (Ill) 
C (graph.) (0002) 
Si (100) 
Si (110) 
Ge(lOO) 
W (100) 




-hco(MeV) 



Figure 8: The parameter k (|72I) versus the photon energy hoj for the crystal channels listed 
in Table [El 



on the distance of a few dechannelling lengths is neghgible. It would be 
substantial if the relative spread SE/E of particle energies would be compa- 
rable to or larger than the ratio X^/Ld. The latter ratio, however, can not 
be made smaller than 10~^ isl fl while modern accelerators usually have a 
much smaller relative energy spread. The same is true for the energy spread 
induced by the stochastic energy losses of the channelling particles due to the 
interaction with the crystal constituents and the radiation of photon. It was 
shown in Ref. 33| that at initial energies of ~ IGeV or smaller, the average 
relative energy losses of a positron in the crystalline undulator AE/E are 
smaller than 10~^. Clearly, the induced energy spread SE/E ^ AE/E is 
safely below the ratio X^/L^. From these reasons, we ignored energy spread 
of the particles in our calculations. 

In conclusion, we have studied the propagation of a modulated positron 
beam in straight and bent planar crystal channel within the diffusion ap- 



■^Note that the corresponding quantity in ordinary ultraviolet and soft x-ray FELs, the 
inverse number of undulator periods l/N^ = Xu/L, is usually of the order of lO^'^-lO^'* 
[3l| . That is why these FELs are so demanding to the small energy spread of the electron 
beam. 
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proach and presented a detailed description of the used formalism. We in- 
troduced a new parameter, the demodulation length, which characterizes 
the penetration depth at which the beam preserves its modulation. It has 
been demonstrated that one can find the crystal channels where the demod- 
ulation length sufficient for producing coherent radiation with the photon 
energy of hundreds of keV. This opens the prospects for creating intense 
monochromatic radiation sources in a frequency range which is unattainable 
for conventional free electron lasers. Developing suitable methods of beam 
modulation would be the next milestone on the way towards this goal. 



Appendix A. Appendix: Solving equation (1661) . 

Appendix A.l. A series expansion at k <^ 1 

At small values of k, the solution of equation fl66|) can be found in the 



form of power series. 

The Laguerre function Li,{^) (which is a special case of the Kummer 
function, L^{^) = M(— z/, 1,^) can be represented as 

oo p -1 n/2 

(Appendix A.l) 

where J„(. . . ) are Bessel functions and the coefficients An are defined by the 
following recurrence relation: 

Aq = I (Appendix A. 2) 

^1 = (Appendix A. 3) 

A2 = ^ (Appendix A.4) 
[nAn^i - (1 + 2u)An-2] . (Appendix A.5) 



Keeping only the leading term in ( Appendix A.l ) (this approximation is 



valid if 3 <^ 1), equation (166|) can be reduced to 

Jo (^y^(l + 2z/(o))(l + z)jo,iv^) = 0, (Appendix A.6) 

where u^^^ is the zero-order approximation to the root of equation (1661) . 
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Equation ( Appendix A. 6 ) is satisfied if 
y^(l + 2z/f )(l + z)jo,iv^ 



J0,k, 



jo,fc is a root of tlie Bessel function: Jo(j, 



0,k) 



(Appendix A. 7) 
0. Here the subscript k = 



1,2,3,... enumerates tlie roots of the Bessel function and the corresponding 
approximate solutions of equation (l66l) . Solving ( Appendix A. 7 ) for 
results into 



1 

2' 



(Appendix A. 8) 



Keeping higher order terms in ( [Appendix A.lj ) and expanding the Bessel 
functions around the the zero-order approximation ( Appendix A.8 ), one ob- 
tains a series expansion of Vk{i^)'- 



i^fc(K) 



l-Uo,fc 1 



1 
2 



1 
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4 Jo,i 



« . Jo,fe 
-Jo,i 



".^0,1" 



(a/k)^ -[{Appendix A. 9) 



720 - - 

(Dots stand for higher order terms with respect to k). For the functions (E 
and (l68l) expansion ( Appendix A. 9 ) takes the form 



/5fc(«:) 



Jo,i (jo% - 17jo% + 54) 



Jlk - 2 



+ . . . 



-K Appendix A. 10) 



(Appendix A. 11) 



Appendix A. 2. Numerical solution 

Expansions ( Appendix A.9 ), ( Appendix A. 10 ) and ( Appendix A.ll ) 
fail at K > 1. Therefore, equation f l66|) has to be solved numerically. In 
our numerical procedure, we use the series representation for 1/^(3): 



oo 

E 

j=0 



(Appendix A. 12) 



Equation was solved by Newton's method. At small k the value found 
from the series expansion (Appendix A.9) was used as initial approxima- 
tion. Then k, was gradually increasing. At each step, the equation was 
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solved and the solution was used as initial approximation for the next step. 
During this procedure, the roots //^(k) were slowly moving in the complex 
plane along the trajectories shown in figure Appendix A. 2 The functions 




Re(v) 



Figure A. 9: The trajectories drawn by the roots ^fc(K) of ([66| in the complex plain at 
varying k. The arrows show the direction of motion of the roots when k increases. 



ak and are plotted in Figures Appendix A. 2 and Appendix A. 2 



Appendix A. 3. Asymptotic behaviour at 1 

As it is seen from Figure Appendix A. 2 some of ^^{k) approaches integer 
real numbers as k — )■ oo. This the case for A; = 1, 3, 4, 6. For these solutions, 
the asymptotic behaviour can be found. 

Let us represent Vkii^) in the form 



i^fc(«:) =nk + Skin), 



(Appendix A. 13) 
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Figure A. 10: The function ak{K) obtained by numerical analysis. 



where = z/fc(oo) is an integer number and the function 6k{K,) goes to zero 
at K — 7- oo. 



Substituting ([Appendix A.13|) into ([Appendix A.12|) and expanding around 



6k = one obtains 



E 

.j=o 



(Appendix A. 14) 



j=0 



Here Pn^ih) is a polynomial of the order whose explicit form will not be 
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Figure A.ll: The function I3k{n) obtained by numerical analysis. 



needed in the following. 

At I3I <^ 1, the infinite sum in ( Appendix A. 14 ) can be approximated by 
an integral and evaluated by Laplace's method: 



E 

j=0 



(Appendix A. 15) 



The polynomial PnS'h) (Appendix A. 14) becomes negligible with respect 
to the exponential at large 3. Similarly, the leading order term dominates the 
first sum in ( [Appendix A.14[ ). The asymptotic expression for Lyi}) takes, 
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therefore, the following form 



(Appendix A. 16) 



Using ( Appendix A. 16] ) and taking into account ( Appendix A. 13 ) one 
obtains the asymptotic expression for the root of equation fl66l) : 



1 



[nkl 



2nfe + l 



exp 



1 + t 



(Appendix A. 17) 

This equivalent to the following asymptotic behaviour of the functions fl67l) 
and (EHl) 



akin 



2nk + 1 

Jo,i 

(Joo)!!!^' 
2"^-i(nfc!)2 
2nk + 1 
2jo,i\/^ 



(Appendix A. 18) 



exp 



sm 





2 2 



(Appendix A. 19) 



2"Knfc!) 



exp 



cos 



Uk 

2 2 



It has to be stressed once more, that not all solutions of equation ( l66l) have 
the above asymptotic behaviour. Among the solutions represented in figures 



Appendix A.2||Appendix A.2[ ([Appendix A.17|)-(|Appendix A.TOj) is valid 



only for k = 1, 3, 4, 6 with Uk = 0, 1, 2, 3, respectively. 
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